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Abstract 

We consider entanglement through permeable interfaces in the c = 1 (1+1)- 
dimensional conformal field theory. We compute the partition functions with the inter- 
faces inserted. By the replica trick, the entanglement entropy is obtained analytically. 
The entropy scales logarithmically with respect to the size of the system, similarly 
to the universal scaling of the ordinary entanglement entropy in (l+l)-dimensional 
conformal field theory. Its coefficient, however, is not constant but controlled by the 
permeability, the dependence on which is expressed through the dilogarithm function. 
The sub-leading term of the entropy counts the winding numbers, showing an analogy 
to the topological entanglement entropy which characterizes the topological order in 
(2+l)-dimensional systems. 
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1 Introduction 

Conformal interfaces provide a natural framework to extend the (l+l)-dimensional 
conformal field theory (CFT) with boundaries. Taking into account the important 
roles played by the boundary CFT in condensed matter physics and string theory, one 
expects that the interface CFT also opens up new directions in such studies. In fact, 
interesting properties have already been found: For example, all the symmetries of the 
rational CFT are generated by a class of interfaces called topological interfaces [1]. 
Topological interfaces also transform one set of D-brancs to another [2,3]. For (po- 
tential) applications and a (partial) list of references, we refer to [1-13] and references 
therein. 

When a system consists of two (or more) sub-systems as in the interface CFT, the 
quantum correlation, i.e., entanglement, between the sub-systems is a useful probe to 
the system. The entanglement entropy is a measure of this entanglement. In (1+1)- 
dimensional systems, the entanglement entropy of the ground state shows a sharp 
contrast between the critical and the non-critical regime [14] and, at the critical point, 
there appears a universal logarithmic scaling with respect to the size of the system 
characterized by the central charge [15,16]. In (2-|-l)-dimensional systems, the leading 
term of the ground state entanglement entropy scales linearly as the boundary size of the 
system, whereas the sub-leading term characterizes the topological order of the systems, 
named as the topological entanglement entropy [17,18]. Entanglement is essential also 
in quantum computation and information. The references on the entanglement entropy 
in (l+l)-dimensional systems with defects inchidc [19-21]. 

In this note, we consider the entanglement entropy in the c = 1 interface CFT. 
A class of conformal interfaces in this theory has been constructed which interpolates 
perfectly transmitting and reflecting interfaces [6] . In addition to the physical relevance 
of the c = 1 CFT, because of the fact that these permeable interfaces are simple but 
possessing structures characterized by some parameters, we expect that they provide 
useful insights into the interface CFT. The entanglement entropy may be a useful probe, 
and its role in the interface CFT is of interest. 

In section 2, we briefly summarize the entanglement entropy with slight generaliza- 
tions in the case of the interface CFT. We also introduce the c—1 permeable interfaces. 
In section 3, we compute the partition functions with the interfaces inserted. In sec- 
tion 4, using this result and the replica trick, we obtain the entanglement entropy 
analytically when the sizes of the two CFT's joining at the interface are equal. The 
entropy scales logarithmically with respect to the size of the system, similarly to the 
universal scaling in the case without interfaces. However, its coefficient is not constant 
but controlled by the permeabihty. The sub-leading term counts the product of the 
winding numbers. This shows an analogy to the topological entanglement entropy in 
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(2+l)-dimensional systems. In deriving the entropy, we adopt two approaches: One is 
based on the BernouUi polynomials and numbers, which provides a general method for 
carrying out the replica trick. Another is a direct evaluation of a sum by an integral 
for the large size of the system. In the course of verifying the equivalence of the re- 
sults from the two approaches, we find that the scaling coefficient is expressed by the 
dilogarithm function. In section 5, we conclude with a brief summary. Some useful 
formulas are collected in Appendix. The entanglement entropy in the c — 1 interface 
CFT has been discussed in [22] in a different setting from ours and in the context of 
the boundary entropy and its holographic dual. 



2 Setup 

2.1 entanglement entropy in interface CFT 

We consider two (l-l-l)-dimensional CFT's defined on a half complex plane Rew > 
and Rew < 0, respectively. The interface between CFTi and CFT2 lies along the 
imaginary axis Rett; = 0. The conformal invariance requires the continuity condition 

Ll - L\ ^Ll- L\ , (1) 

at the interface, where -^^ (a = 1, 2; n e Z) are the left and the right Virasoro 
generators of CFT^. 

We are interested in the entanglement entropy of the ground state. The entropy 
is defined by the von Neumann entropy of the reduced density matrix for the ground 
state pi =Tr2 |0)(0| as 



5 = -T:ipilogpi = -^IVipf 



(2) 

K=l 



where Tr^ stands for the trace over the degrees of freedom in CFT^. The trace of the 
K-ih. power of the reduced density matrix is represented by a partition function on a 
X-sheeted Riemann surface TZk whose branch cut runs along the real axis from w — Q 
to 00 [23]: 

'I^i^f = |^ = ^if(iy/^'^exp[-£ dFwm], (3) 



where represents the fields of CFTi and CFT2. The interface is inserted on each sheet 
of TZk- The normalization factor 1/Z^{1) assures Tripi = 1. From this path-integral 
representation, the entropy is given by 

S^{l-dK)\ogZ{K)\^^^. (4) 
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As in the case without interfaces, how mixed the reduced density matrix is depends on 
the correlation across the interface, and hence the von Neumann entropy measures the 
entanglement of the CFT's through the interface. 

To evaluate the partition function with the interface inserted, we move to z = log w 
plane. Introducing cutoffs at |w| = e and \w\ = L [15, 19], the K-sheeted Riemann 
surface TZk is mapped to a rectangular whose lengths along the real and the imaginary 
axis are (log(L/e), 27tK), repsectively. The interface is mapped to Imz = (2m — l)7r/2 
(m = 1, 2K). In the following, we set e = 1/L for simplicity. 

Here, we impose the periodic boundary condition along the Re 2; direction. When 
the interface is absent, this reduces to the treatment in [15], and we can check that 
our final results actually reproduce those in [15] as a special case. (Our L corresponds 
to E/e in [15], as is clear from our setting.) Also, when the interface is topological 
(see below), our partition function is regarded as a generalization (in the c = 1 case) 
of the generalized twisted partition functions discussed in [5]. In order to carry out 
the analysis independently of the boundary conditions, it is desirable to generalize the 
analysis in [15, 16, 24] based on the conformal invariance to our case with interfaces. 

Now, our partition function is a torus partition function with 2K interfaces inserted. 
Recall that the ground state density matrix is represented by the product of the ground 
state wave function \I'o^'o in deriving the path-integral representation (j3]), and that \l/o 
(^q) gives the path-integral on the lower (upper) half w-plane. One then finds that 
the interfaces at Imz = (2m — l)7r/2 for odd m and those for even m are hermitian 
conjugate with each other. Therefore, we find that 

ZiK) = Tn (/i2)^ /12 ■ ■ ■ {luY , (5) 

where q = e^^'^* with t = tt/{2 logL), we have denoted the interfaces by lu and (/i2)^ 
and each of them appears K times alternately. We have also rescaled the rectangular 
in the z-plane so that it becomes of the standard lengths (27r, 27rt). The subscript "12" 
is put to make explicit the fact that I12 joins CFTi and CFT2 in this order. The above 
expression is manifestly symmetric with respect to CFTi and CFT2, as it should be. 

2.2 c = 1 permeable interfaces 

We now specialize our discussion to the case of the c = 1 permeable interfaces. A 
general way to construct conformal interfaces is to use the folding trick [4]: Let us set 
f = Kew. Then, by flipping the sign of f for f < 0, CFT2 comes to live on the f > 
half plane. In the course, the interface becomes a conformal boundary, which can be 
expressed by a conformal boundary state of the tensor product theory CFTi C?) CFT2. 
In fact, the condition of the conformal invariance of the interface ([T]) becomes that of 
the conformal boundary, Ll^ + L^^ — {L^„ + -^-n.) = 0; since the left and the right movers 
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and the positive and the negative modes in CFT2 are exchanged, respectively, by the 
folding. Conversely, unfolding a conformal boundary state gives a conformal interface. 

In this way, from the boundary states of the tensor product theory of two free 
bosons, the c = 1 permeable interfaces are obtained [6] : 



r± (":/3) 



G 



± (q,/3) 



'i2(fei,fc2)i'^±^ - ^i2(fci,fc2)i'^±) n 

n=l 

00 

{"./3) (a - 



G 



12 {kiMy 



'-^12 (fcl,fc2)' 



°° ^ ( S± Q:?_„Q:?_„-5± a?_„a2 _5± q,]_^q,2 +5± 0,25,2 

e 



^+ ^ e^(^"-*^'^) I A;2iV, fciM) (A;iA^, /tsM] 



Af,M=-oo 

CO 



^- e^(^"-*^'^) I fciM, z^aA^) (/^lA^, fcsMi 

Af,Af=-oo 



(6) 



where 



S 



_i_ , =F cos 20-1- — sin IQj^ 
=F sin 20 cos 20± 





kik2 




sm2e± 



1/2 



(7) 



and 



tan 04 



fc2i?2 

kiRi 



tan 0_ 



2k2RiR2 
ki 



(8) 



a^, a" (a = 1, 2; n G Z) are the modes of the free boson 0^ compactified on a circle 
with radius Ra- They satisfy [a^, a^] = m6m,nS"''' and similar expressions for a,". It is 
understood that (n > 0) implicitly act on G^2 from the left and a^, (n > 0) 

from the right, ka is the winding number for (pa- |^a;'"^a) is the oscillator vacuum for 
(pa with the momentum Ua/Ra and the winding number m^. Its dual is denoted by 
{ua, ma\- Ii2 is obtained from the boundary state with one Dirichlet and one Neumann 
boundary condition. 6^ is the angle between the 0i and the Neumann direction in the 
target space. is obtained by T-dualizing the boundary state for Jj^. In unfolding, 
there is a choice of which CFT is unfolded. Changing this choice gives "anti-interfaces" 

[11], which are equivalent to (/]^)^ with some signs of the parameters flipped. One 
can check explicitly that the continuity condition ([1]) is satisfied by these interfaces. 

The matrices S"^ control interactions between CFTi and CFT2. When 6± are a 
multiple of 7r/2, the two CFT's decouple and the interfaces become totally reflective. 
When 6± are odd multiple of tt/4, the interfaces become totally transmissive. In this 
case, the interface is called topological, since each of the left and the right energy- 
monentum tensor becomes continuous across the intreface, i.e., L,^ = and = L^, 
and hence the interface can be freely deformed. Note that the identity operator is 
included as a special case of the topological interface I^2(ii){'^/^)- 
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3 Partition functions with interfaces inserted 



In this section, we compute the partition function with the interfaces inserted, 
i.e., Z(K) in ([5]) for Jj^. The computation in the following can be regarded as a 
generalization of (part of) that in [11] for the fusion of the interfaces. 

3.1 case of 

Here, we consider the case of J^^. To carry out the computation, it is convenient to 
first focus on a unit of the products of the operators in Z{K), 

J = /+g^g+^i(/+)tg^S+^J, (9) 

and rewrite the quadratic oscillator parts in /{^ and (/i2)^ as 

respectively, where c = cos26'_)_ , s = sin2^_)_. We then linearize the quadratic forms by 
an identity 

which is valid when all Ai,Bi {i = 1,2) are commuting with each other. After the 
linearization, one can explicitly put the creation operators {n > 0) on the left 

of or (G^2)^ the annihilation operators a^, (n > 0) on the right. Further 
pushing the Virasoro generators to the oscillator ground states in (^12)^ using 

e""g^o = g-^oge""?!^ a^]2d commuting the creation and annihilation operators between 
G12 and (^^'2)^ one finds that 

nr ^2^* /• j2 T| ^ ^ 

/ vr^ / vr^ 

n ^ 

X JJ^ g-^i:2^nial_„+{c5„i-sz„2)al_„ . . j~J g- ^ Wn2g"a^+(-s»i)„i+ciS„2)g"aJ, ^ ^]_2) 
n n 

where Znt are the components of ^„ etc., and 

G' = glY^ q'w.M+'w^M-i^i^^]^^ k^M) {k2N, kiM\ , (13) 

N,M 

for kik2 with 



For kik2 = 0, we have different expressions of G' due to the change of the zero-mode 
structure in G^2- We separately discuss this case later. We then take the K-th power 
of J, commute the creation and annihilation operators, and perform the Zn2- and Wni- 
integrals so as to maintain the linearity of the oscillators. Relabeling Zni,Wn2 as Zn, Wn, 
we find that 



Z{K) = Tr, J'' = gf J] g^^^^ 



N,k^ M"r'=fei N,knMl ri~ '° 



(15) 
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K 



fe=i 



iV,M=-oo 
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(16) 



with D„ = (1 - c^?^")-^ and z)t 
Since 



V2i?i sin^^y V cos^^ / ' 



kiRiM\2 



.2Risme+J ' V cose+ y ' ^'^^^ 

the sum over N, M gives a product of the theta function -d^. The remaining P„ are 

evaluated as follows. Introducing a 47^^- vector *y = (Re2;i^\ Im^;!^'*, Re , Im 

•), the exponent in P„ is expressed as — ^ - • where is a x 4K symmetric 

matrix 

( U C ■■■ XJ^ 



M, 



K 



XJ u c 

X] u c 



with C 



a(l2 - a^) 



a(l2 + (J^ 



(18) 



In is the n X n unit matirx, cr^ is a Pauh matrix, and a — —s^q^'^Dn/'i, h = cq^{l — 
q^^)Dn. Performing the Gaussian integrals then gives P„ = D^[d.ei Mk\~^^'^ ■ The 
determinant of Mk here is regarded as a generalization of the circular determinant (see 
e.g. [25]), and obtained similarly: 

K 

detMx = JJdet(l+u;fcC + u;^^t') 



fe=i 



K 



Am2 



(19) 
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where Uk = e^'''''/^, 4 = cos{2nk/K) and = (1/2) [^/l - 2(c2 - s2)g2n + g4n ± _ 
g^")]. To derive the last hne, we have used a formula flA.ip . The above expression 



shows that det Mk and hence Z{K) are actually analytic in K. We also notice that 
Y[ Pn gives rise to a product of the theta function -di. With the help of a formula (lA.2p . 
we finally obtain 

Z(K) ^ ,f |.r-.«3(^m;)*3(?^).-'(2a) n .rV.|2«) . (20) 

for kik2 7^ 0, where rjir) is the Dedekind eta function, and 

vrz/fc = arcsin^|s| sin — j . (21) 

The di part is similar to the oscillator part of the amplitude between D-branes at 
angles. This is naturally understood, once we notice that, in the ii"-sheeted Riemann 
surface TZk-, the array of the interfaces resembles pairs of D-branes at angles. 

3.2 case of 

The case of /{^ is similar. For the oscillator part, it turns out that the integral 
expression for P„ is obtained by replacing in flTBl) c = cos 26'+, s = sin26'+ with 
— cos2^_, sin2^_, and thus the final expression by 6^ i.e., -Ri — > l/2i?i. 

For the zero- mode part, the expression corresponding to f|T7j) is also obtained by 
Ri — > l/2i?i. Therefore, Z{K) in this case is obtained from (!20|) by Ri l/2i?i 
(and hence 0+ 9-). This is expected, since /{^ is constructed from the boundary 
state in which CFTi is T-dualized compared with the boundary state for J^^. 

3.3 special cases 

So far, we have considered the case of kik2 7^ 0. When kik2 = 0, while the analysis 
of the oscillator part remains the same, the zero-mode structure and the product of 
Gf2 and (G^g)^ change. Repeating similar computations, one then finds for that 
the product of ^^s's in ((20]) is replaced by Gf = {itkl/2Rl)^^ {2itklRl) for ki ^ 
0,^2 = and = 9^ {itkl/2Rl)ef {2itklRl) for ki = 0,k2^ 0. When both ki and 
k2 vanish, the original boundary states and hence the interfaces are not well-defined, 
since Cardy's condition is not satisfied. We will not discuss this case. From ([H]), one 
also finds that ^1^2 = implies s = (unless taking the decompactified limit (or its 
T-dual) R12 = 0,00 which is not covered in our setting). Z{K) is then simplified as 
Z{K) = g"^ r]'"^^ {2it) , where O = Oi or O2. The results for /{^ is obtained by 
Ri — > l/2i?i as above. 
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4 Entanglement entropy 



Given the partition functions with the interfaces inserted, we would now hke to 
discuss the entanglement entropy. In the following, we concentrate on the case of J^^, 
since the results for /{^ are easily read off from those for 1^2- We also focus on the case 
with kik2 7^ 0, unless otherwise stated. 

To compute the entropy via (jlj), we need the analytic form of Z{K) in K. A way 
to obtain it is to continue the product in (l20l) with respect to K, and another is to 
use the last expression in f[T^ in terms of pj. We first adopt the former with the help 
of the Bernoulli polynomials and numbers. This provides a rather general method to 
carry out the replica trick. We then use the latter, which is more straightforward. In 
the course of showing the equivalence of the results from the two approaches, we find 
that the entropy is expressed by the dilogarithm function. 

We start with the result of Z{K) in fl20l) . Since the modular parameter t = 
7r/(21ogL) is small for L ^ 1, it is convenient to evaluate it by the modular transfor- 
mation r —^/t- One then finds that (when |s| ^ 0) 

Z{K) = (^+1^1)'" ^-(K^3).mt^^iK)/tf^ ^ 0{e-^/')) , (22) 
where is a positive constant and 

k=l 

with < Uk < I (k = 1, K - 1). Note that (fiK) is of the form 

K-l 



^(K) = J2f{j^) , /(x) = ^arccos^dslsinvra;). (24) 

k=l 



Since f{x) is analytic around x = 0, we expand it as /(x) = Ylm=o fm^"^- A useful fact 
here is that Ylk^"^ expressed by the Bernoulli polynomials bn{x) and numbers 6„ 



as in (1A.4I) . From this and properties of bn{x),bn summarized in Appendix, it follows 
that 



oo „ 

K=i ^-^ m + 1 

m=0 



K=l 



After plugging the explicit form of /(x) and changing the variables as m = 
arcsinh (|s| sinhvrx), we apply the result to (jl]), and obtain 

S = a{\s\)\ogL-\og\kM, (26) 
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up to terms vanishing for L ^ 1, where 

a(|.|) = M _ ^ J\(^^l + {\s\/smhuf -l)du. (27) 

We find that the entropy has a logarithmic scahng with respect to the size of the 
system L, but the coefficient cdsj) is a function of It turns out shortly that cdsl) 
is expressed by the dilogarithm function. The sub-leading term counts the product of 
the winding numbers. This is analogous to the topological entanglement entropy in 
(2-|-l)-dimensional systems characterizing the topological order [17,18]. We also note 
that the entropy is a function of 9+,ki,k2 only, and does not depend on a, P, Ri, R2 
explicitly. 

In special cases, cr(|s|) is simplified. First, let us consider the topological case, 
\s\ = 1. Since the identity is included as a special case, this case should reproduce 
the universal scaling of the ordinary entanglement entropy without interfaces. In fact, 
one finds that a{l) = c/3 = 1/3, which agrees with the result [15,16]. Next, when 
|s| is small, one can show that the second term in fl27|) is log and that 

(t(|s|) —>■ \s\/2. This implies that the leading term decreases as \s\ does, which also 
agrees with the fact that the oscillator part of the two CFT's are decoupling as |s| — * 0. 
The result for small | s | is derived also by directly expanding Uk in fl2Tl) . For general | s | , 
one can check that o"(|s|) monotonically interpolates these two cases. This supports 
an intuition that the entanglement changes according to |s|, since |s| is the strength 
of the interaction between the two CFT's. In [9], certain refiection and transmission 
coefficients are introduced as probes of conformal interfaces. For the c = 1 permeable 
interfaces, they give and s^. Compared with those coefficients, one finds that the 
entanglement entropy (!26|) can probe a little more details of the interfaces. 

When the sizes of the two systems are L and A — L, the entanglement entropy 
without interfaces scales as (c/3) log[(A/7r) sin(7rL/A)] + const. [15], where the sub- 
leading constant term is independent of L [16]. In our case with interfaces, the entropy 
should also be symmetric under the exchange of the two CFT's, and the above scaling 
should be reproduced in a special case. A possible form for L 7^ A/2 satisfying these 
requirements is S* = (7(|s|) log[(A/2) sin(7rL/A)] — log |A;iA;2|- 

Here, some comments for special cases may be in order. When kik2 = and hence 
s = 0, it follows from the result in section 3.3 that the entropy exactly vanishes: S = 0. 
This confirms the fact that the two CFT's decouple in this case. When kik2 7^ 
and hence s 7^ (unless in the decompactified limit), the entropy might appear to be 
negative for small enough s. This, however, is not the case: To obtain (122]) by modular 
transformations, we have used t/s^ <^ 1 for t^s's. Thus, the result in fl26|) is valid when 
the first term is large enough. In fact, since L is the cutoff in our setting and can be 
arbitrarily large independently of other parameters, this condition is always satisfied 
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by taking large enough L. Note that, however small s is, the two CFT's couple through 
the zero- modes if kik2 7^ 0. In order to analyze the case where kik2 7^ and slogL 
is small with actually finite L, one may need to develop a method to compute the 
entanglement entropy for finite systems, e.g., by generalizing the results in [16]. 

One can also derive the entropy by applying the expression of Z{K) = 
g'^'ds'dsq'^^^Yl in terms of pj. Recalling the formula (jl]), we first evaluate 



n=l 

-1 

— F{y)dy + F{0), (28) 

TTt Jo 



K=l 

n=l n=l 



as t — *• 0, where 



F{y) = log \s\ — y + y^l + (|s| / sinhy)^ arcsinh ^sinh?//|s| j , (29) 

and F{0) = 1 + log|s|. We have used PnP~ = jspg^", and assumed that s is not 
vanishing so that F{0) is not divergent. In the case without the derivative Ok, a similar 
approximation by an integral is subtle, since the summand is singular at t = 0. Thus, we 
instead note that H^i \k=i ~ e~'^*/^77~^(2it), which after a modular transformation 
gives Er=ilog^n \k=i - 7r/(12t) + log(2t). Together with ^3^3 ^ \s\/ {2Kt\kik2\) for 
small t, the entropy is obtained as S = cr{\s\) logL — log |/cifc2|, with 

1 2 

~<\s\) = -.+- F{y)dy. (30) 

6 TT^ Jo 

Compared with the previous result fl26l) . a{s) should agree with <j{s). To show this, 
we consider their derivatives: 

, 1 2 f°° dw arcsinh , 2 [°° dz arcsinh 2; 
a'(s) = ^ / , a'(s) = / — , 31 

and 

„ —2 f°° dw (w/s^) arcsinh w —2 f°° dz (s^z) arcsinh 2; 



vr^s Jo ^{l + w^){l + wys^y ' vr2s Jq ^{1 + z^){1 + s^z^ f ' 

(32) 

where we have made changes of variables w = sinhw and z = s~^sinhy. The integral 
for cr"{s) here is performed as 



erf s)" = — ( \ / — „ ^ arcsinh w — arcsinh— 



s — 1/s \ \ w'^ + s"^ s 

logs 
s-l/s 



w=0 

(33) 
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One then finds that — (7r^s/2)cr(s)" is also given by the above, namely a{s)" = a-(s)", 
since the integral representations of s • a" and s ■ a" are related by s It is easy 

to confirm that the integration constants are also the same, e.g., by checking special 
values (7(1) = a{l) = 1/3, a'{l) = a'{l) = 1/4, which verifies a{s) = d{s). 
As a by-product, we find by integrating (t"{s) that (t{s) is expressed as 

^(^) = ^ + J + A f(« + 1) M« + 1) logs + («-!) Li2(l - s) + (s + 1) Li2(-s)l ,(34) 

where Li2(-2) is the dilogarithm function. We summarize some properties of Li2(-2) in 
Appendix. Using them, one can rederive the values of cr(l), cr'(l), and the small-s 
behavior of cr(s). 



5 Summary 

We have obtained the partition functions with the c — 1 permeable interfaces in- 
serted, and the entanglement entropy of the corresponding interface CFT analytically. 
The entropy scales logarithmically with respect to the size of the system, as in the 
case without interfaces [15, 16]. Its coefficient, however, is not a constant but a mono- 
tonic function of |s| controlling the permeability, and is given explicitly in terms of the 
dilogarithm function. The sub-leading term of the entropy counts the product of the 
winding numbers. This is analogous to the topological entanglement entropy, which 
characterizes the topological order in (2 -|- l)-dimensional systems [17, 18]. 

Otii results show that the entanglement entropy is a useful probe to the system, as 
in the case without interfaces. It would be interesting to study how general our findings 
are: For example, does the entropy always contain the topological information of the 
system? Does it always show the scaling as in the case without interfaces? Regarding 
such studies, it would be useful to generalize the analysis based on the conformal 
symmetry [15,16,24] to the case with interfaces. A complication with interfaces is that 
one has to keep track of the shape of interfaces under conformal transformations. It 
would also be interesting to consider implications of the entanglement entropy in the 
context of condensed matter physics and string theory. 
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Appendix 

In the main text, we apply the formulas 



" ^ 2r7r 



Yl^x^ - 2xy cos{e + —) + y^^ = x^" - 2x V cos + y^" , (A.l) 

r=l 

We also use the Bernoulli polynomials (see e.g. [25]) bn{x) (ri = 0, 1, 2, ...) defined by 

( ^^^^ -LTTj 



n=0 



Their derivatives are b'^{x) = nbn^i{x). At a; = 0, they give the Bernoulli numbers 
namely, bn = &n(0). with odd index vanish except for bi, and 6o = 1; = — l/2j ^2 = 
1/6,64 = —1/30,... . One also has 6„ = 6„(1) for n 7^ 1, and &i(l) = 1/2. By the 
Bernoulli polynomials, the sums of powers of natural numbers are expressed as 



n-l 

[m , ^] 

k=l 



+ 1) 5^ A;'" = bm+i{n) - W (n, m = 1, 2, ...) . (A.4) 



The Bernoulli numbers with even index have an integral representation. 



00 ^2n.-l 



62n=4n(-l)"/ ^-^dt. (A.5) 







1-e' 



From this, it follows that 



1 



00 ^^m~l 

n=i 2 .Iq 1 — 



dnbm+i{n) , = 5„,o + ^5^^ + - {-^^) / ^—^,dt . (A.6) 



m + 1 

In section 4, we use the dilogarithm function defined by 



00 



T - / \ sr^ z''' riog(l-U7l , ,, 

a"(s) in (l33l) is integrated by using the above integral representation and 

dz U2{z) = z U2{z) - z - (1 - z) log(l - z) . (A.8) 

12 



From Li2(l) = 7r^/6, Li2(— 1) = — 7r^/12, one can check the values of cr(l) and cr'(l). 
To derive the small-s behavior of cr(s), useful formulas are flA.7p and 



77 



Li2(l - z) = -Li2(z) - logzlog(l - -2) + Y • 



(A.9) 
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